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In the spinfoam framework for quantum gravity, we investigate the conditions to have a physical 
quantum state for the Barrett-Crane model for the 4d quantum gravity path integral. More precisely, 
we look at boundary states stable under evolution with respect to the spinfoam kernel. Hee, we 
look at the simplest case of a single 4-simplex boundary and show that the requirement of working 
with a physical boundary state fixes the width of the semi-classical Gaussian wave-packet for the 
boundary 3d geometry. This is directly relevant to the graviton propagator calculations done in this 
framework, since the Gaussian width enters the numerical factors in front of the graviton correlations 
in the large scale asymptotical limit. Finally, we discuss the application of our computations to the 
Barrett-Crane model beyond the first order (of a single 4-simplex in the bulk) and to the more 
recent EPRL-FK spinfoam model. 
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Introduction 



Spinfoam models provide a background independent framework for a regularized path integral for (loop) quantum 
gravity. They define transition amplitudes between quantum states of geometry through state-sum models which 
can be understood as discrete space-time geometries. In this context, a recent development has been a proposal for 
, ^0 , reconstructing the graviton propagator in this discrete setting from correlations between geometrical observables such 
as the areas of elementary surfaces This allows to extract some semi-classical correlations at large scales, which 
we hope to compare with the perturbative calculations performed in quantum general relativity treated as a quantum 
>> ■ field theory. 

Since the original proposal, they has been a lot of work developing this line of research, mainly focusing on creating 
the new mathematical tools needed for these semi-classical calculations and on using these computations as a criteria 
to select spinfoam models with a correct semi-classical behavior and discriminate between them (see e.g [2]). All these 
developments hint towards the fact that the graviton propagator in spinfoam models lead back to Newton's law for 
gravity at large scales while being regularized at the Planck scale. This behavior has been confirmed by numerical 
simulations in the simplest cases [3|, |j] . 

The main ingredients of these "spinfoam graviton propagator" calculations are: a suitable boundary state peaked 
on a classical 3-geometry, the spinfoam amplitudes for the bulk geometry and relevant observables probing the space- 
time geometry. The main spinfoam models used to define the bulk amplitudes are the Barrett-Crane model, which 
exists in both its Euclidean version O L B and its Lorentzian counterpart [1, H|, and the more recent EPRL-FK models 
5^ , [Tol - [T4T | and their generalizations [Til lq | . 

Considering boundary states, while most of the recent work has focused on building quantum coherent states 
with good semi-classical properties, one important issue is the requirement that the states are physical i.e solve the 
Hamiltonian constraint on the 3d boundary. As we will review and discuss later, this "physical state" criteria can be 
entirely formulated in term of a compatibility equation between the boundary state and the spinfoam bulk amplitude, 
which translates the requirement of stability of the boundary state under the evolution defined by the spinfoam model. 
This question was investigated in the framework of the Ponzano-Regge spinfoam model for 3d quantum gravity, more 
particularly in a toy model where the three-dimensional space-time is triangulated by a single tetrahedron [l7|. In 
this case, the correlations between length observables can be thoroughly studied both analytically and numerically 
(l8| . In this context, the requirement of working with a physical boundary state was shown to determine explicitly 
and uniquely the width of the gaussian wave-packet defining the state This width is relevant in the context of 
the geometrical correlations because it enters the exact numerical factor in front of these correlations. Therefore, if 
we would like to have the exact correlations and not only their scaling properties, we need a definite prediction of 
that width. 

In the present paper, we investigate the consequences of the physical state requirement for the Euclidean Barrett- 
Crane model for the simplest case of a space-time triangulation constructed from a single 4-simplex. In this context, 
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we show that this requirement fixes uniquely the width of the quantum boundary state (in term of the classical 
data) similarly to what happens in the 3d toy model. In the next section, we will introduce and quickly review the 
framework of the spinfoam graviton propagator. Then in section II, we will review the definition of the (Euclidean) 
Barrett-Crane spinfoam model. Section III presents our results about identifying suitable physical boundary states 
for the Barrett-Crane model. We consider both the simplest cases of a decoupled Gaussian ansatz and the more 
interesting case of a coupled Gaussian ansatz involving the Hessian of the Regge action (for a single 4-simplex). 
Finally, we will conclude by a short discussion on the relevance of our physical states under renormalization of the 
Barrett-Crane model and on the possible extension of our results to other spinfoam models. 

Before proceeding to the analysis of the Barrett-Crane spinfoam amplitude for a 4-simplcx, we would like to 
discuss to which extent our notion of boundary states stable under spinfoam evolution can be considered as an actual 
physical state. From the strict perspective of the canonical framework, a physical quantum state should solve the 
Hamiltonian constraint (and actually also solve the vector constraints generating space diffeomorphisms). Thinking 
of this requirement in classical terms, this means that the 3-metric on our boundary should be physical, i.e satisfy 
the constraints, and that we should furthermore consider equivalence classes of boundary 3-metrics under the gauge 
transformation(s) generated by the constraint(s). In our spinfoam framework, considering the 3d case with a single 
tetrahedron in the Ponzano-Regge model, there is a single truly physical state from the canonical point of view, which 
solves all the constraints. This is the completely flat boundary state, which is given by the 6j-symbol. The procedure 
followed in is to look for semi-classical states on the same tetrahedron, which are peaked around a certain classical 
3-metric satisfying the constraints (here defined as the lengths of the tetrahedron's edges) and which are stable under 
evolution defined by the Ponzano-Regge spinfoam amplitude. The requirement of being fixing the tetrahedron's edge 
lengths amounts to a gauge- fixing of the Hamiltonian (and diffeomorphism) constraint (s). This allows us to ide ntify 
a single quantum boundary state for each classically admissible edge- length configurations of the tetrahedron [19j . 
Each of these quantum boundary states describe a flat physical 3d-metric within the tetrahedron. Thus we refer to 
these boundary states as physical in the spinfoam context, although they are not strictly physical states from the 
canonical viewpoint. 

We follow the same logic in the present study of the quantum 4-simplex in the Barrett-Crane model. The basic 
spinfoam amplitude for the 4-simplex implies that the bulk 4-metric is flat, although we will discuss in the conclusion 
the possibility of effective spinfoam amplitudes for the 4-simplex inducing curvature. Thus we expect once again a 
unique physical boundary state for 4-simplex from the strict canonical theory definition. Nevertheless, we will look for 
semi-classical states which are peaked on fixed values of the area of the boundary triangles and which are stable with 
respect to the Barrett-Crane spinfoam amplitude. Fixing the triangle areas amounts again to a partial gauge-fixing of 
the Hamiltonian and diffeomorphism constraints. And we then expect a unique such stable semi-classical boundary 
state for each sets of classically admissible areas. These can be considered as particular gauge-fixings of the unique 
canonical physical flat state. By misuse of language, we shall nevertheless refer, in the specific spinfoam context, to 
our quantum boundary states stable under spinfoam evolution as physical states. 



I. SPINFOAMS AND THE "GRAVITON PROPAGATOR" FRAMEWORK 



The spin foam approach to quantum gravity provides a path integral formulation for the canonical framework of 
loop quantum gravity (LQG). It is based on a reformulation of general relativity of an "almost-topological" gauge 
field theory and the spinfoam path integral is constructed as a discretization of the continuum path integral. Without 
focusing especially on the specific case of quantum gravity, we describe below the generic features of the spinfoam 
framework for an arbitrary gauge group G. As it is well-known, the Hilbert space of quantum geometry states for loop 
quantum gravity is spanned by spin network states \tp) where in tp = (7, ji,i n ): 7 is a graph, ji is a "spin" labeling an 
irreducible representation of the gauge group G associated to the link I of the graph, and i n is associated to the node 
n of 7 and labels intertwiners. We consider a 4d space-time region M. with a 3d boundary E. The spin network state 
\ip) defines the quantum state of geometry of the boundary S, then the spin foam amplitudes define the dynamical 
probability amplitude of that state and are supposed to contain the whole dynamical content of quantum gravity. 
More precisely, the standard ansatz for local spinfoam amplitudes can be presented in the following general form: 

K[i/>] = Yl w(C)^l[A f (j f )l[A e (j f ^ e )l[A v (j f , ie ) (1) 

C\dC=i^ jf,i e f e v 

where the sum is taken over two-complexes C and their "coloring" c = {jf,i e }. The 2-complexes are constrained to 
fit the graph 7^ of the spin network i/> at the boundary. The representations jf are associated to the faces f of C and 
the intertwiners i e to its edges e. They are also constrained to fit the representations ji and intertwiners i n of the 
state tp living on the boundary. Then a spinfoam is defined as a colored two-complex, namely a couple (C, c). Finally, 
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the spinfoam amplitude is made of four types of factors. First, w(C) is a statistical weight that depends only on the 
two-complex C (similar to the symmetry factor for Feynman diagrams in standard quantum field theory). Af are 
weight factors associated to the faces and A e are amplitudes associated to the edges of the 2-complex. These three 
types of factors can be interpreted as measure factors of the discrete path integral defined by K[tp]. All the dynamical 
information is encoded in the vertex amplitude, A v (jf, i e ), which is an amplitude associated with each vertex v of the 
two-complex C and depends on the spins jf and intertwiners i e living on the faces and edges around that vertex. 

The amplitudes are usually assumed to be local, that is they depend only on the coloring of adjacent simplicial 
elements. Thus, Af is a function of the representations located on the face /, A e is a function of the intertwiner 
assigned to e and the representations on the faces containing e, whereas A v depends on the representations on the 
faces and on the intertwiners on the edges containing the vertex v. Generally speaking, the choice of the vertex 
amplitude A v corresponds to the choice of a specific form of the hamiltonian operator in the canonical theory. 

Since a spin network state defines a quantum state for the boundary geometry, spin foams are thus interpreted as 
representing the states of the bulk defining the quantum space-time interpolating between given boundary data. This 
can be clearer when considering an initial spin network tp and a final spin network ip'. Then a spin foam T : ip — > ip' 
defined by the 2-complex C bordered by the supporting graphs of ip and ip' , respectively 7 and 7', represents a 
transition from the spin network state ip = (7, ji,i n ) into ip' = (7' , j'i,i' n ). Nodes and links in the initial spin network 
\ip) evolve into 1-dimensional edges and faces. New links are created and spins are reassigned at vertexes. This defines 
a foam-like structure whose components inherit the spin representations from the initial spin network ip and are at the 
end compatible with the spin representations of the final spin network ip'. The propagator kernel K[ip, tp'} is obtained 
summing over all the spinfoams compatible with the boundary data: 

k[v,v']= E A ^^'\ (2) 

where we have introduced here a more compact notation than in the previous formula |T]): the sum over the compatible 
spin foam J- gathers the sum over the two-complexes C compatible with the graphs 7 and 7' and the sums over the 
G- representations compatible with the representations ji and j[ and over the intertwiners compatible with i n and i' n . 
Af[iJj, ip 1 ] is the spin foam amplitude associated to the spin foam J- : ip — > ip' interpolating between the initial and 
final states. It is given as above as a product of vertex, edge and face amplitudes: A v , A e , Af. This special case 
where the boundary is made of two disconnected pieces, which we can interpret as the initial boundary and the final 
boundary, allows a clearer connection with the canonical framework, whose goal is to define transition amplitudes 
between initial and final states. 

In the present work, we will consider the case where the boundary is connected and the kernel K is then defined as 
a function of only one boundary spin network ip: 

K[ip] = £ Arty] (3) 

■?"|a.F=</> 

where we used the same compact notation as above. 

Before going further, we have to emphasize that the sum over spinfoam J- = (C, c) is not exactly well-defined. 
Usually, for fixed C, the sum over colorings c is well-controlled. On the other hand, controlling the full sum over two- 
complexes is a much more subtle issue. It is nevertheless non-perturbatively defined as quantum field correlations in 
the context of group field theory (see e.g. [II])- In this context, Ajr actually depends on the GFT coupling constant A. 
Indeed the statistical weight w{C) for a two-complex C is given by the symmetry factor of the two-complex (considered 
as a Feynman diagram for the GFT) times a factor \ v where V is the number of vertices of the two-complex. Then 
the sum over 4-geometries in K[ip] can be written as a power series in A; that is: 

00 

K[ip] = £ X v K v [tp] 
v=o 

where Jfy[^] = Y]j?v\ v A^jr[ip] is the sum over spinfoams J- v with V vertices. Assuming that the perturbative 

expansion in A is well-defined and that we are working with a suitable semi-classical boundary state ip, we expect 
that the dominant contributions come from the simplest spinfoams and we consider the lowest order terms in the 
expansion of K[ip] in A as giving the leading order for the path integral and correlations. 

Let us now describe the spinfoam framework for deriving the graviton propagator from correlations between area 
observables. We start by considering a semi-classical spin network functional ^[V"] peaked on a classical 3d metric 
q for the boundary E. Let us remind the reader that tp — (7, ji,i n ) labels basis states of the Hilbert space of spin 
networks. We further require that this boundary state vE , g [?/'] induces a space-time structure in the bulk peaked around 
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the flat Minkowski metric. In particular, this normally fixes the classical boundary q to be the 3-metric induced on 
£ by the Minkowski metric on M. Then we construct correlations between the metric fluctuations for the chosen 
boundary state ^> q : 

W abcd (x,y;q) = K W\ (ip\h ab {x)h cd {y) \i/>)¥ q [ip] = Tr [K h ab (x)h cd (y) ¥,] , (4) 

V- 

where the trace is taken over the Hilbert of spin networks. Here x and y are two points localized on the boundary E. 
In our discrete spin network setting, they are usually determined as nodes of the graph 7 underlying the spin network 
state The metric fluctuations h ab (x)h cd (y) are usually constructed as geometrical quantities. They are two basic 
types of such geometric observables in the discretized geometry setting of spinfoams: the diagonal components of the 
metric tensor can be interpreted as areas and the off-diagonal components as dihedral angles between simplices. 

Finally, this formula defines the 2-point function for the gravitational field in the spinfoam framework. It can be 
considered as the equivalent of the standard 2-point function of the conventional quantum field theory framework, 
which defines the graviton propagator [l[: 

W^ pa (x,y) = (0\T{h^(x)h pa (y)}\0)- (5) 

In the spinfoam setting, ^[V'] is a semi-classical state peaked on both intrinsic and extrinsic geometry This 
boundary data will determine the 4-metric induced in the bulk. Moreover, it should be a physical state, that is a 
solution a the Hamiltonian constraint. This "physical state criteria" can be formulated by the two following conditions: 

The first condition is the normalization of the boundary state, while the second condition translates truthfully the 
requirement to work with a physical state and can be considered as the "Wheeler-deWitt" condition. It is the 
consequences of these conditions on ^B q that we are going to study. 

II. THE BARRETT-CRANE SPINFOAM MODEL 

We now focus on a specific spinfoam model: the Barrett-Crane model of 4d Riemannian Quantum Gravity [j| . 

A. Simple Spin Networks and the Barrett-Crane Amplitudes 

The spin foam amplitude for the Euclidean Barrett-Crane model is given by the evaluation of a spin network ip 
with group Spin(4). Spin(4) is the double cover of SO(4). Using the homomorphism Spin(4) = SU(2)l x SU(2) fl the 
irreducible representations (irreps) of Spin(4) are labelled by two halh-integers: X = corresponding to the 

irreps of the two SU(2) sectors and group elements decompose as the product of two left and right rotations G = gL9n- 
In the Barrett-Crane model, each node n of the graph 7 of ip is associated to a Spin(4) group element G n and each 
link I is labelled by a simple irreps of Spin(4) which corresponds to a couple of half-integers 2/ = the SU(2) 

irreps are the same in the left and right sectors. The evaluation then reads Q 

A[i>)= f HdGnHK^Gj^Gw), (7) 

JSpin(4) n ■ 

where we use the Haar measure dG = dg^dgu. The kernel Kx(G) is the matrix element of G on the SU(2)-invariant 
vector |X, 0) in the X representation with I = a simple representation of Spin(4). Here SU(2) is the diagonal 

rotation group, corresponding to the subgroup of 3d rotations. Expressing the invariant vector in term of left/right 
components 

l Z <°> = ^EC-iy-^J.m^lj.-m)* (8) 



where dj =2^ + 1 is the dimension of the SU(2) representation of spin j. It is straightforward to show that the kernel 
K-x is simply given by the SU(2) character Xj(g), defined as the trace of the group element in the j-representation of 
SU(2). Parameterizing the SU(2) group elements as 

g(0,n)=cos(0)l + ism(0)h-a, 0e[O,n], (9) 
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the characters depend entirely on the class angle 9 (half the rotation angle) and are expressed as: 

Then using the properties of invariance of the Haar measure, it can be proved that the relativistic spin network 
evaluation is actually a 3d object involving only integrals over SU(2) |2lj : 



A 



[ YidgnYl—XnigAgtii))- (11) 

JSU(2) V i d 3- 



This spin network evaluation then defines the vertex amplitude for the Barrett-Crane model. 

The Barrett-Crane model was originally derived for specific two-complexes C, which are the dual 2-skeleton of 
4-dimensional triangulations A4. Starting with the triangulation A4, the dual two-skeleton C is constructed by 
associating to each simplex a point in its interior. An edge of the skeleton connects two points, which corresponds 
to two 4-simplices having a common tetrahedron. Then, a triangle of the triangulation A4 correspond to a face of 
the two-complex C formed by the edges of the dual tetrahedra having this triangle in common. This duality allows 
to derive the spin foam amplitudes from a quantization of the geometry of a 4-simplex. Indeed, the amplitude for a 
single 4-simplex defines the vertex amplitude of the spinfoam model. Arbitrary 2-complexes are constructed by gluing 
4-simplices together and the corresponding spinfoam amplitudes are given by the product of these vertex amplitudes 
corresponding to each 4-simplex. 

Considering a single 4-simplex, the vertex amplitude is constructed as the evaluation of its boundary spin network. 
The boundary graph is constructed as the dual of the boundary of the 4-simplex: each tetrahedron is mapped on 
a node of the graph and each triangle shared by two tetrahedra is mapped to a link between these nodes. This 
gives a graph s with 5 nodes n G {!,••• ,5} and 10 links between them, as represented in fig. [T] Then the vertex 
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FIG. 1: The 4-simplex (or pentahedral) boundary spin network. We label the nodes n = 1, ...5. In the dual picture, they are in 
correspondence with tetrahedra of the boundary triangulation. Two of them are represented. The ten links nm, on the other 
hand, are dual to triangles. Consider for instance the link 45: this is dual to the triangle shared by the tetrahedra 4 and 5. 
Associated with the link 45 is the dihedral angle #45 between the tetrahedra 4 and 5. 

amplitude A[s] is given up to normalization factors by the evaluation of the boundary spin networks and defines the 
Barrett-Crane {10j}-symbol: 

5 10 

17 d 9n II x i< (9~^)9t(i)), (12) 

( 2 ) i=\ 

following equation (fTTj) . The normalization ambiguity leads to an ambiguity in the definition of the Barrett-Crane 
model. It corresponds to arbitrary gluing factors between 4-simplices, or equivalently to ambiguity in defining the 
spinfoam edge amplitude A e . We will come back to this issue in the next section. 




B. The 4-Simplex and the Boundary States 



We now focus on the simplest 4d boundary spin network: a pentahedral that we denote s = {T,j a b,i a ) which is 
the dual of a single 4-simplex. We have introduced the notation (a b) for i linking the nodes a and b. a, b = 1 • • • 5 
label the five nodes of the pentahedral; in the triangulation picture, there are in correspondence with tetrahedra of 
the boundary 4-simplex and links i = (ab) are dual to triangles. Thus, in the triangulation picture, while the five 
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group elements g a £ SU(2) are associated to the 5 tetrahedra of the 4-simplex, the ten representations ji = j ab for 
a =/= b can be seen as attached to the triangles of the 4-simplex. The amplitude Ajr of a spin foam T = (C, c) where 
the two-complex C is bordered by the boundary 4-simplex graph V and its color c is compatible with {j a b, i a } is then 
given (up to normalization factors) by: 



A[s] = n^ J 2 II (MdjJ) 1 '* {10j} (13) 

a<b a<b 



where A e (jf) is the edge amplitude left unspecified here and the {10j}-symbol, the vertex amplitude of the spin foam, 
defined by (fT2|) . Let us point out that A[s] = K[s) in the first order in the group field theory coupling constant A in 
the specific case of the pentahedral boundary spin network s. We now focus on the {10j}-symbol which is defined by 
the following group integral: 



{10j} = / [dg a f 5 n^a^ 1 ) (14) 

• / SU(2) a<b 

The {10j} symbol admits a geometrical interpretation associated to the 4-simplex. The key fact is that it can be 
written as an integral over ten class angles 9 a b G [0,tt] of the group elements ga 1 g b - Qab is the dihedral angle and 
we use the convention 6 aa = 0. In the triangulation picture, it is the angle between the outward normals to the 
tetrahedron a and the tetrahedron b of the 4-simplex. The {10j} symbol then reads, 

{10j}= f dfi[6 ab ]l[xj a M, (15) 

J a<b 

where the measure on the ten class angles dfi[9 ab ] takes into account that the group elements g^ x g b are obviously not 
independent. The measure is in fact simply given by a constraint 20] 

MOab] = Y[de ab am6 ab S{detG ab ) (16) 

a<b 

G is the Gram matrix, a symmetric 5x5 matrix defined as G ab = cos 9 a i,. This constraint 6(G) contains all the 
geometric information and allows to relate the {10j}-symbol to the Regge amplitude (which describes discretized 
gravity) for a Euclidean 4-simplex in the large j asymptotics [2(J Hl[ . The independent variables are the areas of 
the triangles of the 4-simplex representing the triangulation or equivalently the dimension of the representation j a f>, 
dj ab = 2j ab + 1, since A ab = lpdj ab . The geometry of the 4-simplex is fixed by given the ten areas A ab = dj ab and the 
dihedral angles can be considered as functions ab (dj cd ) of the areas. Actually, we scale the area dj ab instead of the 
spins j ab , then the asymptotic behavior on a single 4-simplex at the leading order 1 is given by: 



p 



{10j} ~ P(d jab ) cos ( w S R [d jab ] ) + D(d Jab ) (17) 



where Sn[dj ab ] — ^2 a<b dj ab ab . The function P(dj ab ) is a slowly varying factor, that grows as £~ 9 / 2 when scaling all 
triangle areas dj ab by £. D(dj ab ) is a contribution coming from degenerate configurations of the 4-simplex. This is a 
non-oscillating term which has no geometrical meaning. It was found to scale like 1/C 2 and thus dominate the large 
spin limit of the {10j}-symbol. Moreover, it was shown in a recent paper Q that even after having removing the 
dominating contribution of order 1/C 2 , there are additional non-oscillating contributions still hiding the Regge term. 
However, this sick term is negligible in most computations we are interested in, such as in the computation of dU (see 
[H) or as we will see in the computation of ([5]), because it will not match the boundary data induced by ^ q which 
peaks the asymptotic around the non-degenerate semi-classical configuration. 

On the other hand, the unspecified edge amplitude A e in (p~3|) reflects the normalization ambiguity of the {10j}- 
symbol and different choices leads to different versions of the Barrett-Crane model; we will specify this normalization 
factor latter. 

Since in the Barrett-Crane model, the vertex amplitude A v , defined by (|14[) . and therefore the propagator kernel 
K[s] = A[dj ab ], given by (fT3|) . are independent of the intertwiners, we can rewrite the two conditions ((6]), restricting 



1 The fact to scale the amplitude with the irrep dimension dj ab instead of the spin j ab has no effect on the leading order. 
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tp q to be a physical state, as 2 : 

(dj ab )\ 2 = 1 for the normalized condition, (18) 



Jab 



^K(d jab )i){d jab ) = 1 for the " Wheeler-deWitt" condition, (19) 

jab 

with the propagator kernel a function of the dimensions dj ab only 

\a<b / 

with /i and cr undetermined coefficients. The edge amplitude A e is taken into account in this factor /i (Ila<& dj ab )' 7 ■ 
For more generic spinfoam models such as the EPRL-FK models, we need to take into account the intertwiners in 
the definition of the spinfoam vertex and of the physical states. This is postponed to a further work. 



III. LOOKING FOR PHYSICAL STATES 



In the following, we will work on the Barrett-Crane model in the large spin limit; that is we will consider the 
asymptotic formula of the {lOjj-symbol (JT7J) in the propagator kernel formula (|20l) . We now focus on the issue of 
the semi-classical boundary state. The function ty q should describe the boundary value of the gravitational field on 
the boundary 4-simplex. We thus consider a state peaked on the geometry of a regular 4-simplex: q = (dj 01 O). The 
simplest possibility is to choose a Gaussian peaked on theses values: 

^(d e)(dj ab ) = e~^a<b,a<d a ?d{d3 ab -d jQ )(d 3cd ~d k] )+iY, a<b Qd jab ^l) 

The phase of this semi-classical state determines where the state is peaked in the conjugate variables: is the variable 
conjugate to the spin j and it codes the extrinsic geometry of the boundary. a a ^ d is a given ten by ten matrix. It 
depends on dj in such a way that the relative uncertainties of area and angle on this state become small in the large 
dj limit, namely: 

»#.0, ^¥¥^0 Va<6 (22) 

Assuming that the matrix elements ~ a< ^~j" ' n ^ ne large spin limit with a which does not scale with dj a , the 
fluctuation determined by the gaussian state (|2T|) are of the order: 



i/2-l 
jo 



(tt g |Ad 3 -„Jtt g ) ^ < _ (* q \A9 ab \* q ) ^ d -n/2^ (23) 



2 It is often useful for spinfoam asymptotic calculations to do a Fourier transform and go back to group variables (see e.g. [2 ill for the 
application to the spinfoam graviton propagator). In the present case, this is straightforward to do if we choose a factorizable ansatz, 
i.e. a product state of the form ip(ji, ■■,jio) = Yli=i 4>i{ii)- We further simplify it by taking an equilateral ansatz, i.e <f>i(j) = <t>(j) for 
all triangles i = 1..10. Then following 21], we introduce the Fourier transform of the state: 

/(<?)= X>0')xjG?)- 

3 

Then the two conditions for a physical state simply translates to: 

J dg\f(g)\ 2 = l, 

f [dgm] 5 n fidmg- 1 ) = [ [de mn ] 10 8(dct[cose mn ]) Yl fiBmn) = i. 

a<b m<n 
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which restricts n £]0, 2[. 

In the following, we thus focus on a Gaussian state as boundary state such as the matrix elements of the ten by ten 
matrix are given by 

"cd ~ -r- = with n e]0, 2[ and a, 6 G R (24) 

in the large spin limit. We now study the consequences of the two conditions (|18|) and (|19[) on this boundary state. 
That is the aim is now to determine the consequences on of the requirement that our boundary state is a physical 
state. 

A. Semi-Classical States: the Decoupled Gaussian Ansatz 

We first start by considering an additional ansatz for the boundary state. We take a factorized boundary state: 

* q \jab] = H<l>Uab) (25) 
a<b 

where each 4>(j a b) is peaked around the background value q — {dj 01 0) which corresponds to an equilateral 4-simplex. 
In this simplest case, the ten by ten matrix a a ^ d reduces to a diagonal matrix alioxio- Such a factorized boundary 
state has been used in Q since it is up to now the only setting in which numerical simulations can be performed. 
However, this assumption has not been tested yet. Could such a decoupled gaussian state capture a true physical 
state? Could it satisfy conditions ([T5]| and (|T^|) ? 

Two choices for a factorized boundary state have so far appeared in the literature: 

• the Gaussian state [l], Q , where each factor is given by, 

4>(j) = e- Q(d '~ d ™)V e ^ (26) 

where a is a complex number. 

• The Bessel-based state [ill HH j where each factor is given by: 

T f dj \ _ j I dj \ 

<f> B (j) = i H ' ] g+gfrLgJ cos (dj . e) (27) 

/ /o(%)-/ 2djo+ i(%) 

In the large spin limit regime, we focus only on the Gaussian ansatz. Indeed, in the large spin limit the Bessel part of 
(f2~Tf reduces to the Gaussian in (|2"6]) . Therefore, at the leading order, the only difference between (|2"6"]l and (|2~7T) is in 
the phase: the phase in ([26|i is complex whereas the phase in ([27]) is real. A gaussian state with a real phase would 
lead to the same results as a gaussian state with a complex phase. The interested reader can find details concerning 
the case of a gaussian state with a real phase in appendix |A] and we now tackle the issue of defining a physical state 
coming from a decoupled gaussian state with a complex phase. Therefore, we consider a boundary state of the form, 

*{d j0 ,e)\ja b ] =JfU e- a ^- d ^ 2 e- ied ^ (28) 

a<b 

with Af the normalization constant and a € C. We now want to test this assumption using conditions (IT8l) and (|19l) . 
These conditions lead to the two following equations on Af and a: 

1 = E l*«*,,e)C**)| 9 - T^yio E e-*WE.<.<*--<*,)" (29) 

jab {j ab } 

and, 

1 =E^J*(<^e)(i*6) (I[d jab rP(d jab ) £ e- $ WS.< t ^-^)^E o <J^(^)- %) (^- dj0 , 

jab {jab} a<?> e=±l 

(30) 
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The first equation corresponds to the normalization condition (fT5|) for &(d ,&)\jab]- The second equation is the 
" Wheeler-deWitt" condition for * e)W m the simple case K\j a h] = fi (Y\ a<b dj ab ) a {10j}. Solving them allow to 
determine uniquely J\f and a in term of the coefficients (i and p. The analysis is done in the large spin limit regime. 
In equation (|30p . we have already used the asymptotic formulae of the {10j}-symbol. Moreover, we also replace a in 
(|29|) and (|30j) by its asymptotic expression 

a + ib 
a h 

with a, 6 € R and n restricted to belonging to ]0, 2[ in order that the asymptotic behavior of the relative uncertainties 
of the area and angle on this state is correct. Then solving the two obtained equations requires to distinguish three 
cases with respect to the value of the power n € ]0, 2[. Keeping in mind that the standard choice of scaling for the 
Gaussian width a is n = 1, the final result can be stated as: 

Proposition III.l. The requirement on a factorized Gaussian state (28\) to be a physical state fixes the width of the 
Gaussian a for certain values of the power n: 

1. < n < 1: the width of the Gaussian a is uniquely determined and the coefficient a is restricted. More precisely, 

a e M + , 

a= lBF^ b = Q (32 ) 



1ft „\ - dff > i 



a = 4 U - n ) m - " - 5 



where P is the numerical factor in front of the asymptotics of the {10 j} -symbol. 

2. 1 < n < 2: there is no solution in this case. 

3. the standard case n — 1 : A solution exists if and only if 

1 

. Then the value of a can be determined graphically or numerically in term of the value of /iP. For example, 

a ~ 1 andb~ 1.99 for pP = 1.01 10 5 . 

Before proving this result, let us comment on the third case. This last case should be the most natural since 
the semi-classical state is then peaked on the same way on encoding the extrinsic geometry of the boundary and 
on dj = Aq encoding the intrinsic geometry. However, the calculations do not look particularly more natural or 
transparent in that case. This leads to the conclusion that the choice of a decoupled gaussian state to define a 
physical state might be too simple and should be modified. A new proposition is given in the next subsection. Let us 
now give the proof of the results stated above. 



Proof. In the large spin limit, the summation in equations (|29[) and (|30[) can then be approximated with an integral 
and we can write: 

[ d{dj)e -m*)to-*or 03) 



and 



1 ~ *L_ f Yl d[d jab ](l[ d jab yP(d jab ) J2 e- s ( a )S 0<J (^ r «o) a +'E 8 < k [<i )ol (e9» l -e)-9W(<f Jar <ii„) 2 ] (34) 
(2M) J 



a<b a<b e— ±1 

The first integral (|33|) is just a Gaussian integral, which can been integrated directly: 



given a first relation at the leading order between Af and a . 

To evaluate the second integral (|34|) in the large spin limit, we first notice that the Gaussian implies 



8dj ab = dj ab — dj a « 1 Va < 6, 
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thus we can expand the Regge action around dj : 

a<b a<b a<b ■ 7ab a<b,c<d 3ab Jc<i 



The Schafli identity implies that 

dS R 



9d jab 



\d ]= d j0 - 0ab(dj)\d j= d 3o -6 



since in the equilateral 4-simplex, all the dihedral angles are equal to 0. And we introduce the Hessian, 

AT ab _ d2S R | _ 9d ab | 

Cd ~ dd jab d jcd '^ =d « ~ dd jed ld >= d '° ' [ ' 

then: 

^fe b ]~iod, e + e^H ai> + i Yl iY X t %^eE^ + ^ E N cd 6 d Jab sd Jcd . (38) 

a<b a<Cb,c<d a<b a<fc,c<rf 

We replace £.R[dj a6 ] by this expansion in f[34|) : 

1 * jn a<b d[d jab ](Yl a<b d jab yP(d jab ) [e-^^^M-iE^,.,^^-^^^ 



(39) 



The first exponential is a rapidly oscillating term in dj ab which will vanish when we perform the integration over dj ab so 
we only have to consider the second term. Moreover, at the leading order in Sdj ab we can replace (Ila<6 dj ab ) a P(dj ab ) 
by d)°°P(dj ) in the integral (l34l) . 

At this point, we recall that P{dj ab ) grows as £~ 9 / 2 when scaling all triangle areas dj ab by £. So we can write 
P{dj a b = dj ) ~ -572 in the large dj limit, where P > is a constant. The precise value of P can be computed 
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analytically or determined numerically, but whose its explicit value does not affect the present discussion since it 
can be re-absorbed in the parameter /i. Indeed, all the equations in the asymptotical regime will depend on P only 
through the combination /iP. 

Then, once again we have to integrate a Gaussian integral: 

1 - Kdi °l^T )P J II d[8d jab ]e- s ^^ (40) 

where M is a ten by ten matrix defined by 

M? d = aSf d - iN$ a<b,c<d, (41) 

with 8®% = 1 if the two couples of indices are the same and it vanishes otherwise; = ^ ab \d j= d jo was explicitly 
computed in [l| 3 : 

N ab = _^ (42) 
a jo 



( 


-4 


7/2 


7/2 


7/2 


7/2 


7/2 


7/2 


-9 


-9 


-9 


7/2 


-4 


7/2 


7/2 


7/2 


-9 


-9 


7/2 


7/2 


-9 




7/2 


7/2 


-4 


7/2 


-9 


7/2 


-9 


7/2 


-9 


7/2 




7/2 


7/2 


7/2 


-4 


-9 


-9 


7/2 


-9 


7/2 


7/2 




7/2 


7/2 


-9 


-9 


-4 


7/2 


7/2 


7/2 


7/2 


-9 




7/2 


-9 


7/2 


-9 


7/2 


-4 


7/2 


7/2 


-9 


7/2 




7/2 


-9 


-9 


7/2 


7/2 


7/2 


-4 


-9 


7/2 


7/2 




-9 


7/2 


7/2 


-9 


7/2 


7/2 


-9 


-4 


7/2 


7/2 




-9 


7/2 


-9 


7/2 


7/2 


-9 


7/2 


7/2 


-4 


7/2 


-9 


-9 


7/2 


7/2 


-9 


7/2 


7/2 


7/2 


7/2 


-4 



11 



where N is a ten by ten real symmetric constant matrix with all coefficients independent of dj . In particular, 
— — \J\ Va < b. Therefore, M is a symmetric matrix with all its diagonal coefficients equal to a + i\J^-^— ■ 
At this stage we have to distinguish the three cases mentioned above. Assuming that a is of the form a = , 

JO 

with a, b e K, n G]0, 2[. The three cases with respect to the power n S]0, 2[ are: 

1. < n < 1 corresponds to a > -r^. In this case, we can neglect the terms of order -X- compared to a. 

2. 1 < n < 2 corresponds toa< -r—. In this case, a will be negligible with respect to the terms of the order 



3. n = 1 corresponds to a ~ This case should be the most natural case since it peaks in the same way 

. 

the triangle areas of the 4-simplex around the background value Aq = dj and the dihedral angles around the 
background value 0. 



Let us now separately study the three cases to solve equation (1401) . 

1. The first case < n < 1 is the easiest one. 

Indeed in this case, the matrix M ~ + can be approximated by a ten by ten diagonal matrix: M ~ 

al= (a + ib)l since a ~ 3> -j— and the ten integrals are then decoupled. Thus, we just have to compute a 

Jo j ° 

one-dimensional Gaussian integral: 

u = J d[Sdj]e -^r = (q^ 4r / 2 o)^ (43) 

which gives a second equation for J\f and a. Therefore using equation (|3"5j) . we obtain the following equation on 



a: 



2Y23?(a) 4 
Finally, expressing a under the form ^jfe^, < n < 1), we get that 



1 / 7T _ (/iP) 1 / 5 d (2a-9/10)7r 



JO 



a € R+, 

°" = 3 (I -71 ) and CT > h (45) 
« = ^, 6 = 0. 

In this case, a is real and positive. Furthermore, we get a condition on the normalization factor of the spinfoam 
vertex a > 1/5. 

2. The second case 1 < n < 2 implies that a ~ -C 

"j "jo 

In this case, the matrix M reduces to M — —iN. We have to integrate: f d[6dj ai ]exp(i^8dj ai Nfij6dj ed ). 
Recall that for a m x m symmetric invertible matrix A with signature tj(A) we have: 



[dXi] exp 



XiAijXj 



xl \det(A)\ 



In our case, N is real, symmetric and its signature, which is the difference between the number of positive 
eigenvalues and the number of negative eigenvalues, is equal to —2. Therefore, e la ^ N ^ = —i and our Gaussian 
integral is an imaginary number which is not compatible with the first equation on Af (|35l) . So, we cannot have 
a ~ 2±ih w ith 1 < n < 2. 



3. The third case n — 1 corresponds to a = s±ik_ 

"jo 

Then the determinant of M is a complex number: its real part and its imaginary part are polynomials of degree 
10 and of arguments a and b. We thus obtain a complex equation for Af and a: 



ll(d- \(W<?-9/2) p r /J. \(10<r-9/2) p / 10 

(2A0io / Jl^ |e (2A0 10 V det ( M ) 
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Combining this equation with the first equation (1551) that we already have on Af and 3?(a) = a we get a complex 
equation on a and b: 



K^) 5/2 -*")~ 5>p \/s(lb))= <48) 

with M = -r^-. This equation implies a condition on the normalization factor of the spinfoam vertex: 

1 

a = — . 

5 

We can then solve numerically the previous equation by plotting a 3d graph representing the square of the norm 
of the complex number given by the left-hand side of the previous equation (|48l) in terms of a and b using Maple 
(see graph fig. [2]). For example, the value of this norm vanishes for 

a~ land b ~ 1.99 for fiP = 1.01 10 5 . 

Therefore, there exists a specific value a for which ^(d jo ,&) is a physical state. However, whereas this case for 
which ^(d jg ,&) is peaked in the same way around the intrinsic and extrinsic geometry (dj ,Q) should be the 
most natural, it is quite complicated. 

□ □ 




FIG. 2: We plot the norm squared of the left hand side of eqn (|48[) for fiP = 1.01 10 5 in term of a and b. The second graph is 
a slice of the 3d graph for fixed a = 1 and variable b. 

The alternative is then to work with a more complication semi-classical state. Indeed although a factorized gaussian 
wave-packet is up to now the only one which has allowed to perform numerical simulations, the previous analysis seems 
to show that it is too simple to catch all the features of a physical state. We show in the next section that a tensorial 
Gaussian state is more adapted to describe a physical state. Indeed we will see that such a state allows to compensate 
the imaginary part which comes from the second derivative of the Regge action and given by the matrix iN and 
consequently to simplify the resolution of this third case studied above. This is what's commonly expected from the 
more general perspective of Regge calculus [24| ■ 

B. The Coupled Gaussian Ansatz 

To go beyond the decoupled Gaussian ansatz, our new assumption is to consider a boundary state of the form: 

ip q [d jab ] = i e -E«<i c <^°5(' l ^r' i J o)( i Jd-' 1 Jo) e 'eE a < l 'i lai (49) 

where a is now a ten by ten complex matrix. Furthermore we choose 

< b d = /3(d h )C + ^ (50) 
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where /3 € 1 and the imaginary part of a is now the conjugate variables of the dihedral angles of the tetrahedron in 
the semi-classical regime introduced in ([37| : 

AT ab _ dOab | 
Cd ~ ddiJ di=d >°- 

Jed 

N£d is the Hessian of Sr and we will see that it is this choice which allows to simplify the construction of a physical 
semi-classical state peaked in the same way on the extrinsic and extrinsic geometry of the 3d boundary. We recall 
that N depends on d,- n such that N = with N a matrix with constant coefficients. 

Proposition III. 2. For < n < 1, we take as ansatz for the width 8 = -jpr (at least at leading order for large jo) of 

the Gaussian state {49\j . This defines the real part of the matrix a through (|50|) . Then the width is uniquely defined 
and the coefficient a is restricted. More precisely, 

=^r g ' (5i) 



a = ^ (j? — n) and a > 



5 > 



Therefore, the case n = 1 appears now in the continuity of the case < n < 1. This can be considered as an 
improvement compared to the previous case of a factorized Gaussian semi-classical state. However, we have not 
solved the case 1 < n < 2, which admitted no solution in the factorized Gaussian state. The resolution of this case 
needs the knowledge of the next-to-leading order of the asymptotic expansion of the {10j}-symbol, which hasn't yet 
been computed. 

Proof. Regarding the asymptotic behavior of 8, the real part of the matrix a, we distinguish the three same cases as 
in the decoupled gaussian state case. 

1. The result of the first case, B ~ -%r with < n < 1, is not modified, therefore, B ^> -j— and the matrix M 

a j a 30 

can be approximated by a ten by ten diagonal matrix M ~ B I and the state is again factorized and we obtain 
the same result as in (|3"2"|) : 



B G K+ 



(? 

a = 3 (I - n) and a > | 



(52) 



The second case is when 1 <n< 2 in /) ~ (a € R). Since B <^ j— in this case, the real part of a 

is negligible compare to its imaginary part. Consequently since this imaginary part compensates the second 
derivative of the Regge action in the exponential of the second condition p^|) . we should go to the next order in 
8dj. To compute this expansion in Sdj we need to know the next-to- leading order of the asymptotic expansion 
of the {10j}-symbol. We can thus say nothing now in this case for the moment. 

The third case, 8 ~ (a G K), is greatly simplified compared to the equivalent case for the decoupled 
gaussian wave-packet. Indeed, applying condition (TT5|) and condition (flT))) to this state, we obtain the two 
following equations for Af and a: 



1 * ^ (j d[6 dj ]e- 2 ^ )2 y = (^) 5 

,10<r-9/2 D / „ ,2\ 10 jIOct— 9/2 „ 



2 lr W 



(a lu ^2\ iu jIOct— 9/2 „ . . 



(53) 



where we use to compute the second equation the same analysis as in the factorized gaussian state case. The 
key element which simplifies everything is that the second derivative of the Regge action added to the imaginary 
part of a is null since a imaginary part is the Hessian. This system of equations can then be simplified in a 
single equation for a: 



1 /^oV_ 4o 0ff " 9) WV ^o NM ' 

2 10 I 2a ) 2 20 



(54) 
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which implies a condition on the normalization factor of the {10j}-symbol in the spinfoam vertex: 



and the unique solution for a: 



a=\ (55) 



a= . (56) 



Therefore, we now get the same value of a in this case as in the first case. Moreover, a is now uniquely determined 
since n = 1. 

□ 

This analysis on the consequences of conditions (fT5|) and (fT9"| on a coupled gaussian wave-packet (|4"!)|) is therefore 
more convicting on the ability of such a state to capture a true physical rather than a too simple state defined as 
a factorized gaussian state. Indeed, the width a of the coupled gaussian state is now uniquely determined by the 
requirement of a normalized physical state, as expected from the 3d case [l9j . 



Conclusion and Outlook 



The analysis done shows that it is possible to determine a physical semi-classical state for the Barrett-Crane model 
in the case of the simplest triangulation given by a 4-simplex. The requirement of a normalized physical state 
determines uniquely the Gaussian width (in term of the precise normalization of the {10j} vertex amplitude). This 
analysis shows also that we cannot take a too simple semi-classical state. For e.g. we have seen that a decoupled 
gaussian state which is peaked on the geometry of an equilateral 4-simplex does not seem to be able to capture a 
true physical state. A coupled gaussian state appears as defining more naturally a physical semi-classical state. We 

recall that in all our calculations, we have used K\[s] the lowest order term in the expansion of K [s] in the group field 
theory coupling constant A and that s symbolizes the 4-simplex boundary graph. The total propagator kernel K[s] is 
given by K[s] = Y^v=i ^ V Ky[s]. This operator should describe a unitary "evolution". In the spin foam framework 
the notion of evolution is not well-defined, however we can require the following normalization condition on K[s] 

K[s]K[s] = (K[s]f = 1 (57) 

that describes the process of creation of a 4-simplex starting from a null 4-volume following by the annihilation of 
this 4-simplex into a null 4-volume again. This condition also says that K is by definition a physical state. Using the 
A expansion of K[s] = X)v=i -^^VHj we can expand this normalization condition in power of A 

oo 

^A y K Vl K V2 =l (58) 

v=i V 1 ,V 2 /V 1 +V 2 =V+l 

which simplifies at the leading order in A in 

\ 2 K 1 [s] 2 = l (59) 

meaning that approximatively XKi [s] is a physical state and introducing a new constraint on XKi [s] = 
•V* Ila<6 dj ab {10j}. It is this new constraint which allows to determine //: 



1 

(u a<b {I0j} 2 



(V) 2 = — , , 2ffr , (60) 



where a > 1/5 in order that \x is finite. This restriction on the domain of validity of a is consistent with the results 
obtained previously. Since up to now we have only considered the lowest order term in the expansion of K[s] in A, we 
have in fact fixed the simplest bulk triangulation and taken into account only a finite number of degrees of freedom. 
The next step would be, keeping the same boundary spin network s, to explicitly write the next terms: K2, K^,... of 
the A-expansion of K [s] = J2v=i A y ify[s]. This is not obvious. A priori, K4 could be determined from K\ by doing 
a 5-1 move. Otherwise, from an effective field theory point of view we can write K [s] = Ki[s] + k[s] where k takes 
into account all the possible counter-terms and therefore all the possible bulk geometries for a given boundary spin 
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network. This is equivalent to look at a coarse-grained lattice in which all vertices can be considered as shrunk to 
a single effective vertex. k[s] is then the weight associated to this effective vertex. In this weight k[s], we expect to 
get a term proportional to the original Barrett-Crane spinfoam vertex amplitude, e{10j}, which would come from the 
contribution of the flat 4-simplex. Moreover, we expect additional contributions representing curved 4-simplex. These 
extra terms can be written in a condensed way as p J dG P(G) Tioj(G), where rioj(G) is the spin network functional, 
on the 4-simplex boundary, evaluated on the non-trivial holonomies G. The flat 4-simplex is then simply the term 
rioj(I) = {10j}, and P{G) is the measure factor giving the relative weight of every curved configuration. Other terms 
will certainly appear in a more precise analysis. The issue is then to determine which terms of k[s] contribute in 
the requirement of a physical state for ip q ; that is, we should understand how the condition J^. K[dj ai ]ip(dj ah ) — 1 
is explicitly modified by the term k[dj ab ] of K[dj ab ]. We would expect that the requirement for the boundary state 
to be a physical state selects the bulk triangulation. From this perspective, the choice of the phase seems crucial: 
indeed it is the choice of the flat dihedral angles in the phase factor of the boundary state that is used to select the 
flat 4-simplex at leading order. The choice of a different phase for if? q should select another bulk triangulation T, 
which might correspond a non-trivial space-time curvature. A detailed investigation is required to confirm this rough 
analysis and to really understand the dominant terms in the computation of the two-point function of the gravitational 
field. Up to now, we have discussed the case of the simplest 4d triangulation defined by a single 4-simplex. We should 
go beyond and study more general backgrounds defined on more complex triangulations. We would then expect a 
more complicated dependence on the boundary spins of both the Gaussian matrix and the phase factor, as happens 
in Regge calculus (see e.g. [24j). To deal with such more complex backgrounds, we would however need to develop 
more powerful analytical tools to probe the spinfoam amplitudes and their asymptotics. We nevertheless hope to 
tackle this issue soon. Finally, another step would be to deal with the EPRL-FK models instead of the Barrett-Crane 
model to define the propagator kernel K[s]. Since the EPRL-FK vertex has also for asymptotic the exponential of the 
Regge action, we are expected to get similar results in the restricted case of a single 4-simplex. However, the behavior 
under renormalization should be quite different since the space of intertwincrs in the EPRL-FK model introduces new 
degrees of freedom in the 3d space geometry compared to the Barrett-Crane geometry which is describe by a unique 
intcrtwiner and seems to have too many frozen degrees of freedom. These degrees of freedom will have to be taken 
into account when gluing 4-simplices together. 



Appendix A: Physical States with a Real Phase 



We will discuss here the case already mentioned of the Bessel-based factorized boundary state: 

i>q[jab] = Y[ ^Btiab) 

where 



a<b 



T t djg \ _ T f d jp \ 

<t>B{3) = . - COS (dj9) 

V/o(%)-/ 2djo+ i(%) 

We are interested in the large spin limit regime and its dj — > oo limit behaves as a Gaussian peaked around d 



(Al) 



(A2) 



Jo- 



d jo TT 



1/4 



exp 



2d 



-{dj-d JO ) 2 



Jo 



cos(djO) 



(A3) 



The difference with the case studied previously of a Gaussian peaked around dj is the phase which is real here. We 
recall that the factorized boundary state assumption has been made in order to perform numerical simulations and 
the choice of a real phase has been done in the work concerning the area correlator to turn it into an exact group 
integral [2l| and perform exact analytical computations. 

We therefore consider now a factorized Gaussian state with a real phase: 



^q[dj*t]j{ IX exp(— a(cf Jai , - d ]0 f) cos(d J0 6) 



(A4) 



a<b 



to see if the conditions necessary to have a physical state are modified. Conditions (|18|) and (|19|) give the two following 
equations for Af and a. 



Jab 



1 



-2Sf(a)E a<(j (^ ai) -^o) 2 



n 

a<b 



1 + cos(2d iab 6) 



(A5) 
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and, 



Jab 



5( ti io (l -' 1 jo) 2 + i X;„<6K aD ( £e ^+»)e)-9(a)(d ja 



{jab} 0<6 



e=±l;?7=±l 



_(A6) 

where we have already used the asymptotic formulae of the {10j}-symbol. In the large spin limit, the summation in 
the two previous equations can then be approximated with an integral and we can write: 



1 



1 



27V 2 



d(dj)e- m ^^- d ^ 2 



1 1 



,i2di@ 



1 



and 
1 - 



(2/V) 



' J a<b a<b e=±l;i7=± 



,-SR(a) J2 a<b (d jab -d n ) 2 + l J2 a<b [d Jab (e8 ab +r i e)-Q( a )(d jab -d :jn ) 2 



(A7) 



(A8) 



The first integral (IA7[) after some changes of variables is just three Gaussian integrals, which can been integrated 
directly: 



N 2 



2K(a) 



1 +e 2 «(=) cos(<2 io 6) 



(A9) 



but we recall that a oc ^r- with n e]0, 2[ from equations (|23[) and therefore in the large spin limit the second term of 

the right-hand side of this relation vanishes and we obtain the same relation at the leading order between M and a 
as in the case of a Gaussian state with an imaginary phase (up to a factor 1/2). That is: 



N 2 



231(a) 



(A10) 



To evaluate the second integral (|A8[) in the large spin limit we expand as previously the Regge action around dj 
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a<b,c<d 



where 
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ing term in dj ab 


will vanish when preforming the integration 


i (IA8II which are 


two Gaussian integral: 



(d . o) (lQ g -9/2)p 

(27V) 10 



HdlSdj 



a<b 



Sd 3ab M? d 8d 3cd + e Sd Jab Qf d Sd jad 



(A12) 



where M and Q are both ten by ten matrices defined by: = a5f d - iN™ and Qf d = a5f d + iN™ with a < 6, c < d 

J cd 



and — 1 if the two couples of indices are the same and it vanishes otherwise. Therefore, M is a symmetric matrix 



with all its diagonal coefficients equal to a + i\ t-t— and Q is a symmetric matrix with all its diagonal coefficients 
equal to a — i\/T -r— . Then the analysis is the same as the one done in the case of a Gaussian state with an imaginary 



phase; we have three case to consider: 
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1. a -r—, that is we consider here that a is proportional to d- n with < n < 1. We have then that M ~ a/d 
and Q — a/d and once again we write a = then we get the same result as in the case of the imaginary 

JO 

fael 

phase (up to some factors 2 in a) (see equations 152")) . that is: < a = j (| — n) and <r > i 

[ a = 32P 4 vr 

2. a <C that is we consider that a is proportional to d~ n with 1 < n < 2. We have to integrate 

"jo • /0 

J d[<5dj ab ] exp(j Sdj ab N^Sdj cd ) + exp(~iJ2Sdj ab N^Sdj cd ). The signature of the matrix N is equal to —2; 
therefore, e* ^^* = — i = — e J<T ^ ^7 and then the previous integral is null so in this case we also cannot have 
a — -£r with 1 < n < 2. 

JO 

3. a ~ -T-: it is the most natural case which peaks in the same way the triangle areas of the 4-simplex around 

a,j 

the background value Aq — dj and the dihedral angles around the background value 0. But once again due 
to the fact we have chosen a factorized boundary state, this case is very complicated although the imaginary 
parts of the determinant of M and of the determinant of N will compensate. This confirms that a factorized 
boundary state is not the most natural to capture a physical state. And with a real phase it is not even possible 
to consider a tensorial state which could compensate iN and —iN like we did in the case of the imaginary phase. 
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